Small Angle Shubnikov-de Haas Measurements in Silicon MOSFET's: the Effect of 

Strong In-Plane Magnetic Field. 
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Measurements in magnetic fields applied at small angles relative to the electron plane in silicon 
MOSFETs indicate a factor of two increase of the frequency of Shubnikov-de Haas oscillations at 
H > Hsat- This signals the onset of full spin polarization above Hsat, the parallel field above which 
the resistivity saturates to a constant value. For H < Hsat, the phase of the second harmonic of 
the oscillations relative to the first is consistent with scattering events that depend on the overlap 
instead of the sum of the spin-up and spin-down densities of states. 



PACS numbers: 71.30.+h, 72.20.My, 73.40.Hm, 73.40. 



A great deal of interest has recently been focussed on 
the anomalous behavior of two-dimensional (2D) systems 
of electrons and holes whose resistivities un- 

expectedly decrease with decreasing temperature, behav- 
ior that is generally associated with metals rather than 
insulators One of the most intriguing characteris- 
tics of these systems is their enormous response to mag- 
netic fields applied in the plane of the electrons |@,||j9) 
or holes (|Jl^: the resistivity increases sharply by more 
than an order of magnitude, saturating to a constant 
plateau value above a magnetic field H^at- 

In this paper we report studies of the resistivity of 
the 2D electron system in silicon MOSFETs in magnetic 
fields applied at small angles (f> with respect to the plane. 
This allows a study of the Shubnikov-de Haas (SdH) oscil- 
lations in perpendicular fields sufficiently small that the 
orbital motion has a negligible effect on the response to 
the in-plane component of the magnetic field. At small 
tilt angles (j>, the SdH oscillations plotted versus filling 
factor have twice the period below Hsat compared with 
the period above Hsat- This implies that the electron 
system is fully spin-polarized in high fields, H > Hsat, 
where the resistivity has reached saturation. Detailed 
examination of the oscillations in fields below Hsat indi- 
cates unusual behavior consistent with electron scatter- 
ing that depends on the product rather than the sum of 
the spin- up and spin-down densities of states. 

Two sihcon MOSFETs with mobilities fi « 
20,000 V/cm^ s at T = 4.2 K were used in these stud- 
ies. Contact resistances were minimized by using samples 
with a split-gate geometry, which permit high densities 
in the vicinity of the contacts while allowing independent 
control of the density of the 2D system under investiga- 
tion. Standard AC four-probe techniques were used at 
3 Hz to measure the resistance in the linear regime us- 
ing currents typically below 5 nA. Data were taken on 
samples mounted on a rotating platform in a ^He Ox- 



ford Heliox system at temperatures down to 0.235 K in 
magnetic fields up to 12 T. 
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MAGNETIC FIELD (Testa) 

FIG. 1. Resistivity of the 2D electron system in silicon MOSFETs 
versus in-plane magnetic field for different densities, as labeled. 
Data are shown at T = 245 mK. The arrow indicates Hsat for 
electron density Us = 0.80 X 10^^ cm . 



Measurements were first taken with the plane of the 
sample oriented parallel to the magnetic field [|ll|. The 
resistance, Rxx, is shown in Fig. 1 as a function of field 
for different fixed gate voltages spanning densities be- 
tween ris = 0.8 X 10" cm-2 and ris = 2.88 x 10^^ cm'^ 
(the zero field critical density for the metal-insulator 
transition is ~ 0.84 x 10^^ cm"'^). Consistent with 
earlier findings ||7|,p|, p^|p^ jl^ , the in-plane ((/> = 0) mag- 
netoresistance rises dramatically with increasing field and 
saturates above a density-dependent field H satins) p^ ]. 

The sample was then rotated to make a small angle 4> 
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with respect to the field, so that the in-plane component 
was almost equal to the total field H\\ k, H, while the 
projection in the perpendicular direction, H± « (jjH, re- 
mained relatively small even in high fields. R^x and R^y 
were measured simultaneously as a function of magnetic 
field for fixed angle </>, temperature T, and density Us. 
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Filling Factor v 
FIG. 2. Resistivity R 

XX versus electron filling factor f — 
7ls$0 / H± showing SdH oscillations for different electron den- 
sities Tig, as labeled. Data are shown at T = 0.25-?^ for two 
different angles (j) between the magnetic field and the electron 
plane: (a) (j) = 6°; (b) (j) = 7.8°. 

For various different densities Ug, Figs. 2 (a) and 2 (b) 
show the resistance Rxx vs filling factor ly at two differ- 
ent angles (f> between the magnetic field and the 2D plane. 
Similar curves were obtained at other small angles. The 
filling factor v = ns^Q/H± was calculated using the re- 
lation Us = H±/ {Rxyec). The Hall resistivity Rxy and 
Hall coefficient Rh were determined from the low-field 
data, Rxy = Rh x H±, {i. e. in fields below the onset of 
quantum oscillations). 

For the lowest densities shown in Fig. 2 (a) at angle 
= 6°, the system is in the high- field saturated regime 
above Hgat for filling factor i/ < 10. The large arrow 
indicates Hgat for Ug — 1-54 x 10^^ cm~^. Quantum os- 
cillations are clearly evident superimposed on the large 
plateau value of Rxx at small v. In this region the period 
of the SdH oscillations corresponds to a change in filling 



factor Aiygat — 2 (including the two-fold valley degen- 
eracy for a 2D layer of (100) silicon). For higher filling 
factors > 10 where the 2D electron system is below sat- 
uration {H < Hsat ), the period of the SdH oscillations is 
twice as long, namely Av — 4. Similar behavior is shown 
for a bigger angle in Fig. 2 (b), where the larger perpen- 
dicular component gives rise to stronger SdH oscillations; 
here the period doubling is found above ^ « 8. 

The period Ai' = 2 of the oscillations a,t H > Hgat cor- 
responds to complete spin polarization of the electrons. 
There are several scenarios that could account for a tran- 
sition to full polarization in strong magnetic fields, de- 
pending on the nature of the ground state of the system 
at = an issue that is currently under debate. 
Here we restrict the discussion to a simple model within 
a single- particle description. Using this approach we 
were able to explain the doubling of the frequency of the 
SdH oscillations at > Hgat- However, the detailed be- 
havior of the oscillations in small perpendicular fields is 
not fully consistent with this model. 

As shown schematically in the inset to Fig. 2, the 
spin-up and spin-down electron bands are split by the 
Zeeman energy Az — gfisH, while the spacing between 
the Landau levels, hujc = heH±/mc, is determined by 
the perpendicular component of the field. We consider 
the progression as electrons are added to the system: for 
small densities, Ep < Az, electrons are added to Landau 
levels in the spin-up band only, corresponding to a SdH 
periodicity Aiy = 2 (including a factor 2 for the valley 
degeneracy in silicon); at high densities, Ep > Az, twice 
as many electrons are required to fill both spin-up and 
spin-down Landau levels, yielding the double period, or 
Au = 4 1 15|. An equivalent argument holds for fixed den- 



sity as one reduces the magnetic field. Thus, the shorter 
period Ai' = 2 at > Hsat signals the onset of full po- 
larization of the electron system |l^. The relationship 
between Hsat and complete spin polarization was also 
found by Okamoto et al using a different experimen- 
tal method. 

The observed period-doubling is consistent with this 
simple model only if the spin-up and spin-down levels 
are degenerate or nearly so, so that a — AzjfiLo^ = i 
with i an integer. The double period should revert to 
a single period when a ~ i + 1/2, corresponding to a 
spin-up Landau level between two spin-down Landau lev- 
els. The ratio a can be varied experimentally by chang- 
ing the angle 4>, or by using the fact that the electron- 
electron interaction-enhanced g-factor (and thus Az) de- 
creases with increasing electron density in silicon MOS- 
FETs ||l^,|l^ . By taking data over a broad range of den- 
sities, we were able to smoothly vary a by more than 
1. Close examination of the data shows that the double 
period in fields below Hgat, although stable over a broad 
range, does break down in a narrow range of densities 
that is different for different angles (/>, as expected within 
this model. 
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SdH oscillations reflect changes in electron scattering 
due to periodic oscillations of the density of states at the 
Fermi level |l8|. For the weak perpendicular fields used 
in our expriments, there is strong scattering and the SdH 
oscillations are small |l9| . Unlike the situation that pre- 
vails in high magnetic fields, where the Landau levels 
are sharp and well-defined, the density of states is best 
represented in this regime by a harmonic expansion [[l7[ : 



Di ^(E) = Dail + ex cos{2Tr{E ± Az /2)/hujc) + Oie^)) 

(1) 

Here E, Az = gi-isH, and huc are the energy, Zeeman 
energy and cyclotron energy, respectively. The small pa- 
rameter e = 2exp{—'K / (lOct)) <C 1 |17| is proportional 
to the Dingle factor Small variations in the re- 

sistivity are proportional to small variations in scatter- 
ing: Ap/p — AW/W. Using the Born approximation, 
W J 5{E - EF)D{E)dE, one can show that: 



e X cos{t:Az /hujc)cos{2TTE'^/hujc), Az < 2El 



cosiiTTE^p/hiVc), Az >2E'l 

(2) 

Here the Fermi energy Ep is measured from the bot- 
tom of the band at iJ = 0, Ep is the Fermi energy 
a.t H ~ 0, and D{E) is the total density of states: 
D{E) = D-f{E) + Di{E). We assumed T = and ne- 
glected higher harmonic terms of order in Eq.(l) as 
well as higher order corrections due to oscillations of the 
Fermi energy. This demonstrates that the SdH period 
changes by a factor of 2 when Az > 2Ep, corresponding 
to full polarization of the electrons. 

The term A ~ ecos{TTAz/fiu!c) of Eq. (2) depends on 
the ratio a — Az/hujc, which is fixed for a given angle 
and electron density and does not vary with magnetic 
field. It determines the overall amplitude of the oscil- 
lations at fields below Hsat, when Az < 2Ep. This 
amplitude has a maximum when a — Az/hujc = i is 
an integer, corresponding to spin up and spin down den- 
sities of states oscillating in phase, and vanishes when 
a = « -I- 1/2 (see Eq.l). Figure 3 shows data over a nar- 
row region near a = z -I- 1/2 where the amplitude of the 
first harmonic is small, allowing detailed examination of 
higher harmonic terms (see Eq.(l)). 

Based on the usual assumption that the SdH oscil- 
lations are determined by the total density of states, 
D{E) = D-^{E) + Di{E), one expects and generally ob- 
serves |20) the progression illustrated schematically in 
Fig. 4: a minimum which becomes progressively deeper 
develops at the center of each maximum (see curves (a) 
and (b)), gradually splitting it into two separate maxima 
(curve (c)). Thus, the minima of the second harmonic 
(curve (c)) are at the positions of the maxima of the first 
harmonic (curve (a) in Fig. 4). 
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FIG. 2. In magnetic fields below saturation, H < Hgat> the 
SdH oscillations are shown for a narrow range of electron den- 
sities near a = Az/fuVc = gpsEl /TlLOc = {i -\- 1/2); due 
to the e — e enhancement of the g-factor, each density corre- 
sponds to a slightly different value of a. Note that the maxima 
of the second harmonic (Ug = 4.14 X 10^^ cm^'^) are in phase 

with maxima of the first harmonic(71s — 3.83 and 4.51 X 10^^ 
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FIG. 2. Schematic of SdH oscillations for a density of states 
D(E) = D-^(^E) + D^[E) when: (a) the spin-up and spin- 
down Landau levels are degenerate, a = Az/fll^c = (t>) 
a = i + 6 for small S; and (c) a — i -\- 1/2. Note that the 
maxima in (a) and (c) are 180" out of phase. Curve (d) illus- 
trates maxima for the second harmonic in phase with maxima 
of the first obtained from assuming a scattering probability of 
the form W = f{D^{E) X Di{E)), in agreement with the 
experimental behavior shown in Fig. 3. 
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However, careful examination of the data of Fig. 3 
shows that the behavior below Hsat observed experimen- 
tally in silicon MOSFETs is quite different: no minima 
develop within the maxima, splitting them into two; in- 
stead, the maxima simply diminish in amplitude and new 
neighboring maxima appear and grow in amplitude. The 
maxima of the first [us — 3.83 and 4.51 x 10^^ cm^^) 
and second [ug — 4.14 x 10^"'^ cm^^) harmonics arc in 
phase (as in curves (a) and (d) of Fig. 4) rather than 
180° out of phase. The origin of this unusual behav- 
ior is not clear and warrants further careful study. In- 
terestingly, the phase relation between first and second 
harmonics observed in our experiments can be obtained 
within the single particle model used earlier if one con- 
siders the product of spin-up and spin-down densities of 
states, W = f{D^{E) x Di{E)), rather than their sum; 
curve (d) of Fig. 4 is the result of such a calculation. 
This suggests there is a sizable contribution to the elec- 
tron scattering from events that depend on the overlap 
of spin-up and spin-down densities of states, perhaps re- 
flecting enhanced scattering of electrons of opposite spin. 

In summary, measurements of small-angle Shubnikov- 
de Haas oscillations indicate that the period of the oscil- 
lations changes by a factor of two at the magnetic field 
Hsat above which the resistance has reached saturation. 
We attribute the abrupt change in period to the onset of 
full polarization of the electron spins. The period dou- 
bling in fields below Hgat is stable with respect to the 
angle between the magnetic field and the 2D plane, and 
is observed for all electron densities except in a narrow 
interval, where the amplitude of the first harmonic of the 
SdH oscillations vanishes and the second harmonic is ob- 
servable. The phase observed for the second harmonic 
relative to the first is consistent with SdH oscillations 
due to scattering events that depend on the overlap in- 
stead of the sum of the spin-up and spin-down densities 
of states at the Fermi level. This unusual behavior may 
reflect the importance of many-body interactions in the 
2D system. 
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